Nature Ecology & Evolution
Preprints posted in the last 7 days, ranked by how well they match Nature Ecology & Evolution's content profile, based on 131 papers previously published here. The average preprint has a 0.10% match score for this journal, so anything above that is already an above-average fit.
Zambelli, F.; D'Addese, G.; Marti-Baena, Q.; Sardanyes, J.; Aguade-Gorgorio, G.; Sole, R.
Show abstract
Genomic instability is a major driver of tumor evolution, promoting diversification and adaptation while simultaneously increasing the accumulation of deleterious alterations. How tumor populations balance these opposing effects remains poorly understood. Here, we introduce a computational framework that explicitly represents diploid genomes, functional gene classes, point mutations, and chromosome-segregation errors in spatially constrained and well-mixed tumor populations. We identify a viability boundary separating sustained tumor expansion from instability-induced population collapse. Within the viable regime, mutation and selection generate a stable distribution of genomic-instability classes that is accurately captured by an analytical replicator--mutator description. Near the viability boundary, tumor dynamics exhibit prolonged extinction transients and strong sensitivity to stochastic fluctuations, with important differences between solid and liquid architectures. Chromosomal alterations further modify growth by creating transient benefits through increased gene dosage and genetic redundancy, while ultimately increasing genomic fragility. Finally, simulated interventions show that eliminating low-instability subpopulations or increasing the global mutational burden can displace tumors beyond their viability boundary and trigger irreversible collapse. These results identify genome instability as both an evolutionary advantage and an intrinsic vulnerability, providing a quantitative framework for developing therapies that exploit the limits of tumor evolution.
Lopez-Idiaquez, D.; Satarkar, D.; Sheldon, B. C.
Show abstract
Most evidence of the consequences of climate change in natural systems has focussed on shifts in mean temperature (1,2), but the effects of extreme climatic events (ECEs) remain far less understood. This is particularly true for very severe ECEs that may occur only once every few decades. Understanding the consequences of these severe events for natural populations is nonetheless critical, since their frequency is predicted to rise under current climate change (3). Here we combine a unique long-term dataset spanning almost five decades of breeding (>20,000 events) and morphological data (>120,000 observations) in adult and nestling great tits (Parus major) and blue tits (Cyanistes caeruleus) with fine-scale temperature records to examine the effects of an unprecedented heatwave in May 2026 on breeding success and morphology. Average temperature during the heatwave (22-29 May 2026) was 7.85 C above the historical record, reaching +10.5 C (+4.32 SD) at its peak (25-26 May). These record-breaking temperatures significantly reduced adult breeding success and nestling bmass relative to expectation in the absence of a heat-wave. Given the heatwave was widespread (Fig. 1A), our findings from a single, exceptionally well-studied population are likely to generalise to other species exposed to the same event, providing key evidence that severe ECEs can substantially harm wild populations.
Kuntner, C.; Philippe, C.; Vraka, C.; Zachhuber, L.; Wanek, T.; Friske, J.; Weissenboeck, V.; Helbich, T.; Hacker, M.; Tanaka, E.; Otsuki, L.
Show abstract
Tissue regeneration requires careful allocation of metabolic resources, yet how organisms adjust this allocation in response to varying amounts of tissue loss remains poorly understood. Here, we show that the regenerative metabolic response is not fixed: the size of an injury regulates how glucose is allocated at both local and organism-wide levels. We first demonstrate that tail regeneration requires glucose metabolism in the axolotl (Ambystoma mexicanum), a salamander capable of regenerating centimetre-scale tissues. We then mapped glucose uptake in axolotls regenerating from small or large tail injuries using positron emission tomography/magnetic resonance imaging (PET/MRI) and the radiolabelled glucose analogue [18F]FDG. Glucose uptake was elevated in regenerating tails compared to uninjured tails. During early regeneration, larger injuries induced higher glucose uptake than smaller injuries, correlating with faster regenerative outgrowth. Larger injuries also increased glucose uptake in distant organs, indicating a systemic metabolic response. Together, our findings suggest that metabolic responses tuned to injury size underlie faithful tissue regeneration and establish PET/MRI as a powerful approach for studying whole-body metabolic dynamics in large regenerating vertebrates.
Clark, A. G.; Jiang, J. Y.; Chitale, M. D.; Cosgrove, E.; Van Elgort, A.; Jain, A. M.; Kelso, J. C.; Cui, X.; Yapici, N.; Lin, C.-c.
Show abstract
Prions, once mainly studied for their pathogenic roles, are now gaining recognition as adaptive elements in microbial physiology. Over one-third of wild yeast isolates harbor prion proteins, yet their impact on host-microbe interactions remains poorly characterized. Given the ecological dominance of yeasts in the Drosophila mycobiome, we leveraged the Drosophila melanogaster-Saccharomyces cerevisiae system to investigate how the mycobiome-derived prion, [MRPL10+], modulates host physiology. We show that flies exposed to [MRPL10+] yeast exhibit significantly enhanced cold tolerance and increased locomotor activity. This effect persists with heat-killed yeast and diluted culture, suggesting a stable, potent bioactive factor. Using the genetically diverse Drosophila Global Diversity Lines (GDL), we identified natural variation in responsiveness to [MRPL10+] yeast. Genome-wide association and functional RNAi screening revealed a gut-brain signaling axis involving genes critical for digestion, intercellular communication, transcription regulation, and neural transmission. Notably, serotonin and octopamine pathways were essential for [MRPL10+]-induced changes in cold tolerance and locomotion, implicating neuromodulatory circuits in prion-mediated microbial signaling. Our findings establish a mechanistic link between a fungal prion and host metabolic and neural adaptation. This work provides the first genetic dissection of a prion-mediated host-microbe interaction, laying the groundwork for investigating beneficial prions in complex microbial communities and highlighting a new dimension of the mycobiomes influence on animal physiology.
Tassios, E.; Pyrgelis, N.; Rinker, D.; Tzermpou, E. M.; Hittinger, C. T.; Rokas, A.; Nikolaou, C.; Vakirlis, N.
Show abstract
Genes encoding novel protein sequences are a ubiquitous feature of genomes. They fuel molecular and cellular evolutionary innovations and frequently contribute to species-specific characteristics. We are now unravelling the processes by which they originate, including de novo from noncoding sequences and through extreme divergence, yet how much and what types of novel proteins evolve through each process is still unclear Does the mechanism of origination shape the structural and functional potential of the resulting proteins? Here, we conducted a broad computational investigation of genetic and protein novelty at the scale of the entire subphylum of Saccharomycotina yeasts. We detected more than 5,000 robust de novo genes across 332 species and compared them to more than 10,000 novel genes resulting from extreme sequence divergence, revealing two distinct modes of evolution of novelty. A remarkable 40% of de novo proteins are predicted to localize to mitochondria compared to only 15% of divergent, with the latter also being substantially longer and more disordered. A detailed analysis of conservatively predicted tertiary structures of novel proteins shows that "invention" of novel folds can happen through both processes but is more likely to occur de novo. We also illustrate cases of evolutionary "re-invention" of existing protein folds from non-coding sequences. Our work deepens our understanding of the origins and importance of novel proteins opening new directions for further structural and functional characterization.
Aga, O.; Moyo, S.; Ferno, J.; Manyahi, J.; Kibwana, U.; Löhr, I.; Langeland, N.; Blomberg, B.; Johnston, I.
Show abstract
Infections with bacteria which have evolved multi-drug resistance (MDR) cause millions of deaths worldwide. Large-scale efforts are gathering genotypic and phenotypic data on MDR bacteria, but methods for learning the structure, diversity, and predictors of evolutionary pathways to MDR have yet to take full advantage of these data. Here, we use evolutionary accumulation modelling (EvAM), an emerging class of machine learning methods with roots in cancer progression, to infer these evolutionary pathways across ESKAPEE pathogens (seven bacterial species that dominate health burdens), using a database of over 635k genotyped phenotypic observations from around the world. We identify global patterns in MDR evolutionary pathways, remarkably shared across multiple ESKAPEE species. Species-specific deviations from these stereotypical pathways are connected with geographical and demographic covariates, facilitating predictions of future MDR evolution. We verify these predictions with several hundred new phenotypes from ESKAPEE samples spanning decades of clinical infections in sub-Saharan Africa, demonstrating the capacity to forecast future MDR evolution from these inferred shared pathways.
Jilani, A.; Allgeyer, E. S.; Li, X.; Guo, M.; Sevilgen, D. S.; Ball, A.; Xiong, F.; McLaren, S. B. P.
Show abstract
The symbiosis with photosynthetic dinoflagellate algae enables corals to build and sustain reef ecosystems. Individual coral polyps hold algal symbionts in their epithelial endoderm cells and lose them under environmental stress, leading to coral bleaching. How the host integrates symbionts into its body plan is not well understood. Here, using a combination of high-resolution imaging, quantitative analysis, and environmental perturbations in the sea anemone Exaiptasia diaphana (Aiptasia) and reef-building coral Pocillopora damicornis, we uncover a spatial organisation of symbionts along the aboral-oral axis of cnidarian polyps that emerges under the long-range translocation of symbionts between host cells through a fluid-filled cavity. The symbiont distribution becomes specifically enriched in the tentacle bud endoderm during Aiptasia polyp morphogenesis. This pattern can form in darkness and with algae-sized inert spheres, suggesting an innate host-intrinsic mechanism. Symbiont-occupied host cells are mechanically constrained within the endoderm and thus unable to rearrange; instead, they go through cycles of symbiont expulsion and re-uptake via the host gastric cavity, with regionally biased rates of these behaviours providing a route to enrich symbionts in the tentacles. Symbiont organisation is remodelled under increased light in adult coral polyps, with a characteristic pattern of reduced tentacle enrichment, lateral clustering and retention in the body column emerging over a timescale of days. Together, our findings reveal that the spatial organisation of symbionts is dynamically regulated in cnidarian host tissues, a capacity that may shape both the establishment of symbiosis and its resilience under environmental change.
Dutta, S.
Show abstract
Bulk telomere length measured from tumour sequencing is routinely interpreted as a property of the cancer cells. However, a tumour specimen is a mixture, and the patient who supplies it has a telomere length of their own. Here I re-analyse published pan-cancer telomere estimates and ask how much of a tumour's telomere length is patient-specific. A calibration step comes first. Whole-genome and low-pass estimates recover the known cross-sectional attrition of leukocyte telomeres with age, at 26.6 bp per year in blood normals, whereas whole-exome estimates do not. After adjustment for cancer type, sequencing centre and sex, the exome slope is minus 0.6 bp per year. In 684 blood-normal aliquots sequenced by both assays, the whole-genome estimate declines at 38.9 bp per year, whereas the exome estimate from the same DNA shows no detectable decline. The difference between assays is 41.5 bp per year, with P = 3 x 10^-10. Because exome data constitute 78.6% of the original resource, downstream analyses use only whole-genome and low-pass libraries. Within those data, tumour telomere length tracks the patient's matched-normal telomere length. The Spearman correlation is 0.395 in TCGA, with positive associations in 22 of 23 cancer types. This finding replicates in PCAWG using a different telomere estimator, with a correlation of 0.472 and positive associations in all 24 histologies examined. Adjustment for cancer type, sequencing centre and library type leaves a regression coefficient of 0.385. The association is also stable after adjustment for age, sex, tumour purity, leukocyte fraction, ploidy, sequencing coverage and continental ancestry, with coefficients ranging from 0.406 to 0.429. Pure normal-cell admixture is rejected as the sole explanation. Under a two-compartment mixture model, the coefficient for host telomere length is expected to equal 1 and the host-by-purity interaction to equal minus 1. These restrictions are jointly rejected with P = 0.001. Tumour purity, leukocyte fraction and age each explain only about 1 to 3% of within-cohort variance and do not alter the cross-cancer ranking. By contrast, the between-cohort coefficient is not directly interpretable. Its apparent near one-to-one relationship with tissue-associated telomere length depends strongly on which tissue supplies the matched-normal reference and on the statistical spread of that predictor, falling to 0.44 when organ-matched solid tissue is used. Bulk tumour telomere length is therefore a composite phenotype containing a replicated patient-specific component. Telomere biomarker studies should include matched-normal telomere length as a covariate rather than treating tumour telomere length as exclusively tumour-intrinsic.
Vigna, A.; Harrouard, J.; Miot-Sertier, C.; Loegler, V.; Marullo, P.; Friedrich, A.; Schacherer, J.; Peltier, E.; Albertin, W.
Show abstract
Brettanomyces bruxellensis is a yeast species associated with diverse fermentation environments and characterized by extensive genetic diversity, including diploid, autotriploid, and allotriploid lineages resulting from independent hybridization events. These lineages are associated with distinct ecological niches and provide a framework for studying metabolic trait evolution in complex genomes. Nitrate assimilation is a relatively uncommon trait among yeasts and has been reported in B. bruxellensis, but its distribution and evolutionary history within the species remain poorly understood. Here, we combined phenotypic characterization of 151 strains with genomic analyses of 946 whole-genome sequences to investigate nitrate assimilation. Growth assays revealed that nitrate assimilation is widespread but unevenly distributed across genetic lineages, with some populations largely retaining the trait whereas others have frequently lost it. Genomic analyses identified extensive variation affecting the nitrate assimilation gene cluster composed of YNR1, YNI1, and YNT1. Nitrate assimilation was strongly associated with both gene copy number and predicted gene functionality, with nitrate-assimilating strains generally carrying more functional copies of the cluster. Leveraging the complex genomic architecture of the species, we independently analyzed primary and acquired genomes in allotriploid lineages and uncovered contrasting evolutionary trajectories following hybridization. While nitrate assimilation genes were generally maintained in primary genomes, acquired genomes showed a higher prevalence of gene loss and predicted loss-of-function variants, revealing asymmetric dynamics between subgenomes. Altogether, our results suggest that nitrate assimilation represents an ancestral trait that has been differentially maintained across B. bruxellensis lineages through a combination of copy number variation, gene degeneration, and genome-specific evolutionary dynamics. These findings provide new insights into how genome architecture and polyploid evolution shape the maintenance and loss of metabolic traits in an industrially relevant yeast species.
Klein, C. A.; Koerkel-Qu, H.; Raya, E.; Guzvic, M.; Irlbeck, C.; Mederer, T.; Spitzl, D.; Czyz, Z.; Schunicht, L.; Seitz, S.; Roth, J.; Rack, B.; Harbeck, N.; Kurdieh, H.; Mayr, R.; Burger, M.; Robold, T.; Hofmann, H.-S.; Weber, M.; Maak, M.; Janssen, K.-P.; Huecker, S.; Kirsch, S.; Werner-Klein, M.; Perry, A. C.
Show abstract
Disseminated cancer cells (DCC) in non-metastatic carcinoma patient bone marrow (BM) are predictive of metastasis. Those detected by epithelial cytokeratin or EpCAM expression have poorly-characterized transcription profiles due to their extreme rarity: 1~2 cells per two million BM cells in every third non-metastatic patient. We here characterize the transcriptomes of DCCs. Single-cell RNA-sequencing (scRNA-seq) of 864 EpCAM-positive cells (from 1,151 cancer patients) in BM or lymph nodes (LN) revealed plasma, immune, myeloid, erythroid progenitor cells and two candidate DCC populations, termed M0-DCC and M1-DCC. M0-DCC, mostly from non-metastatic M0-stage patients, displayed the highest known adult stemness scores, and were transcriptomically reminiscent of human cleavage-stage, preimplantation embryos. M1-DCC represented cancer cells undergoing the epithelial-mesenchymal transition (EMT), corresponding to later, implanting and gastrulating embryos. Detection of early-embryo-like DCC categorised patients at highest risk for metastatic progression. Furthermore, high M0-DCC scores predicted the metastatic potential of human cell lines from the Cancer Cell Line Encyclopedia. M0-DCC gene expression profiles can be reversibly induced from M1-DCC-like cells in vitro. The close correspondence between gene expression profiles in immediate early embryonic development and metastatic founder cell candidates provides strong evidence that the onset of cancer and metastasis recruits mechanisms employed in fertilization.
Perina, F. J.; Thomas, V.; Ketehouli, T.; Mudiyanselage, S.; Jain, M.; Schlathoelter, I.; Goss, E.; Martins, S. J.
Show abstract
Antibiotic-induced disruption of plant-associated microbiomes has the potential to alter host health beyond the directly exposed generation, yet whether the effects of dysbiosis are transmitted through the seed microbiome remains unknown. Here, we investigated the intergenerational impacts of streptomycin-induced dysbiosis in tomato (Solanum lycopersicum), demonstrated that seed microbiome transfer (SMT) restores progeny microbiome function and disease resistance, and characterized the underlying physiological and genetic mechanisms. Parental streptomycin exposure altered the composition of progeny rhizosphere bacterial communities, reduced expression of defense-associated genes, and increased susceptibility to Xanthomonas perforans. Suppression of immune gene expression was strongly associated with increased disease severity, indicating that parental dysbiosis impaired progeny plants ability to mount effective immune responses. Transfer of the seed microbiome from healthy plant donors partially restored rhizosphere community composition, reduced disease severity and recovered defense gene expression of three genes. Together, our findings demonstrated that antibiotic exposure microbiome disturbance generates intergenerational legacy effects that influence plant immunity and disease susceptibility and seed microbiome transfer can counteract this dysbiosis across generations.
Mezawa, Y.; Kumegawa, K.; Morita, K.; Yang, L.; Hirakuri, K.; Yamashita, K.; Shirakihara, T.; Sasaki, R.; Onagi, H.; Kutomi, G.; Maruyama, R.; Orimo, A.
Show abstract
Tumor-promoting myofibroblastic carcinoma-associated fibroblasts (myCAFs) are induced by activation of transforming growth factor-{beta} signaling. However, the molecular basis of myCAF-specific transcriptional programs regulated by TGF-{beta} signaling remains poorly understood. Using a meta-analysis of single-cell RNA-seq data from 132 human breast tumor and non-tumor tissues, we show that myCAFs activate gene regulatory programs relevant to skeletal and cardiovascular development that are associated with poorer outcomes in breast cancer patients. Of note, distal-less homeobox 5 (DLX5), a master transcription factor for skeletal development, is activated in human breast myCAFs at both epigenetic and transcriptional levels. DLX5 expression is also initiated by TGF-{beta}1 treatment in human mammary fibroblasts. Immunoprecipitation and CUT&RUN assays using DLX5-expressing fibroblasts demonstrate that DLX5 interacts with Smad2/3/4 proteins, enabling their cooperative occupancy at shared genomic binding sites of target genes, thereby promoting canonical TGF-{beta} signaling and the myCAF state. DLX5-primed myCAFs also enhance paracrine TGF-{beta} signaling and neuropilin-2 expression to promote collective tumor invasion. Our findings indicate that DLX5 induces myCAF formation and promotes breast tumor progression in collaboration with canonical TGF-{beta} signaling.
Cauldron, N. C.; Dort, E. N.; Weeks, G.; Rogers, D.; Cuomo, C. A. A.
Show abstract
Drug resistance emerges repeatedly in outbreaks of Candida fungal pathogens, but little is known about its origins or persistence. Here, we investigated the evolutionary processes shaping echinocandin resistance in Candida auris, a globally emerging and predominantly clonal fungal pathogen. Genome-wide association across over 600 isolates identified mutations in the {beta}-1,3-glucan synthase gene FKS1 as the most significant driver of resistance to an echinocandin drug. Ancestral reconstruction of this population traced shared resistance mutations among small groups typically consisting of 2-3 closely related isolates, but clusters could include up to 16 isolates. Nearly all resistant clusters consisted of isolates collected in the same year and region, consistent with local transmission. To further examine population-level selection, we measured adaptive signatures in FKS1 and the highly diverged paralog FKS2 across 22,000 genomes. This revealed excess nonsynonymous polymorphisms in FKS1, primarily due to independent, recurrent mutations at resistance hotspots, consistent with parallel evolution and incomplete fixation of adaptive alleles. In FKS2, there is no evidence of hotspots and little support for diversifying selection. Together, these results indicate that resistance mutations emerge under strong genetic constraint, with adaptation restricted to only one FKS homolog and predominantly at mutational hotspots.
Baruah, G.; KC, Y. K.
Show abstract
The shape of density-dependence governs species persistence, and ecosystem stability. Yet, whether per-capita growth declines sublinearily, or superlinearily with density remains hotly debated. Growth rates across the tree of life have been shown to decline sublinearly with density, whereas theory founded on resource competition predicts the opposite. Here, we resolve this discrepancy and show that sublinearity can readily emerge from geometric constraints on consumer interactions. By linking inter individual spacing, movement and interference rates, we derive two limiting-interference regimes, one of which the well-mixed limit recovers the form of classic Beddington DeAngelis interference response. We then developed an individual-based model from first principles which reproduces the derived sublinearity response, and further use empirical data from published consumer-resource experiments that also bears the signature of sublinear density-dependence. Further, embedding the interference mechanisms underlying the emergence of sublinear density-dependence in coexistence theory opens a new regime for species coexistence where classical theory fails to predict. Our framework indicates that non-consumptive interactions are not merely a correction to resource competition but might be a distinct axis along which diverse communities may potentially coexist.
Dupre, G.; Pouget, B.; Martinez-Pineda, A.; Foret-Lucas, C.; Bessiere, P.; Chretien, D.; Ducatez, M.; Vialaneix, N.; Hoede, C.; Marquet, R.; Gaspin, C.; Volmer, R.
Show abstract
High pathogenicity avian influenza viruses (HPAIVs) emerge from H5 and H7 low-pathogenicity avian influenza virus progenitors through mutations that introduce a multibasic cleavage site in haemagglutinin. Although nucleotide insertions recurrently generate this motif, the molecular determinants of insertion and whether particular HA sequences are genetically predisposed to evolve toward HPAIV remain unknown. Combining experimental virology and thermodynamic modelling, we show that insertions arise through polymerase slippage controlled by local product-template duplex thermodynamics within the viral polymerase catalytic site. Predicted RNA secondary structures outside the polymerase are not required for high-frequency insertions and only modestly modulate insertion rates. We formalize this mechanism in HPAIVpredict, which predicts insertion profiles, recapitulates intermediates associated with documented HPAIV emergence events and identifies H5 and H7 sequence backgrounds predisposed to acquire functional multibasic cleavage sites.
Zerbato, B.; Taverna, G.; La Chimia, M.; Pontoriero, M.; Lombardi, S.; Taglietti, L.; Deng, K.; Perrone, G. C.; Hakkola, S.; Vuori, A.; Syriala, T.; De Billy, E.; Barabino, S. M.; Bragato, C.; Pierri, C. L.; La Ferla, B.; Urbanucci, A.; Scumaci, D.; Chiaradonna, F.
Show abstract
Pancreatic ductal adenocarcinoma (PDAC) exhibits profound metabolic rewiring and strong resistance to DNA-damaging therapies, yet how metabolic pathways regulate genome maintenance remains poorly understood. The hexosamine biosynthetic pathway (HBP) integrates nutrient availability with protein glycosylation through production of UDP-GlcNAc, but its role in DNA damage response (DDR) regulation is unclear. Here we show that inhibition of the HBP enzyme phosphoglucomutase-3 (PGM3) reduces DNA repair capacity in pancreatic cancer cells. Transcriptomic and functional analyses reveal that the selective PGM3 inhibitor FR054 amplifies gemcitabine-induced replication stress, disrupts ATR-CHK1 and ATM-CHK2 checkpoint signaling, and selectively impairs homologous recombination. Glycoproteomic profiling identifies the AAA+ ATPase RUVBL2 as a key metabolic-DDR node. Gemcitabine increases RUVBL2 O-GlcNAcylation, with Thr81 identified as a modified residue within the Walker A nucleotide-binding motif. Structural modelling predicts that Thr81 O-GlcNAcylation stabilizes the RUVBL1-RUVBL2 complex without compromising ATP-Mg engagement. PGM3 inhibition and Thr81 mutation similarly reduced ATR and ATM abundance and promoted persistent DNA damage, supporting a role for RUVBL2 Thr81 O-GlcNAcylation in sustaining checkpoint signalling and genome stability. Consequently, PGM3 inhibition induces a BRCAness-like state that sensitizes pancreatic cancer cells to PARP inhibition, both in vitro and in vivo, as well as to ionizing radiation. These findings reveal a nutrient-sensitive mechanism linking protein glycosylation to genome maintenance and identify HBP-dependent DNA repair as a potentially actionable vulnerability in pancreatic cancer.
White, J. R.; Robinson, J. D.; Doremus, M. R.
Show abstract
Heritable bacterial symbionts are pervasive in terrestrial arthropods, often imposing reproductive manipulations to promote their own spread within host populations. Co-infections are common, potentially allowing symbiont co-infectors to hitchhike through a host population. However, adverse thermal conditions can disrupt these communities, particularly when co-infectors vary in their thermal sensitivity. We used a multi-generation experiment to test whether warm (29 {degrees}C) conditions disrupted spread of heritable symbionts through uninfected populations of the spider, Mermessus fradeorum. We tested two common infection combinations: a single infection with a cytoplasmic incompatibility (CI) inducing Rickettsiella or a feminizing co-infection that included a feminizing Wolbachia, the same Rickettsiella, and up to three apparent hitchhikers (two additional Wolbachia strains and Tisiphia). We initiated replicate populations with 1/3 of one infection type and 2/3 uninfected spiders, evaluating population infection rate over 5 spider generations under different temperature regimes. Under cool (21{degrees}C) conditions, Wolbachia feminization drove co-infection to 88% and Rickettsiella CI drove single infection to 83% of host populations. Vertical transmission for all symbionts was high (97-99%) and hitchhiking symbionts also spread effectively. Under warm conditions, feminization and CI efficacy were reduced, and symbionts suffered variably reduced vertical transmission. Warm conditions ultimately destroyed the co-infecting symbiont consortium and impeded symbiont spread. On its own, though, Rickettsiella was still able to increase, despite reduced strength of CI. We hypothesize that contrasting tensions between feminizing spread of the symbiont consortium versus environmentally driven loss of function and transmission may explain observed patterns of mixed infections in field populations of this spider.
Kavanagh, D.; Steel, A.; King, H. E.; Vieira, H. G. S.; Kumar, K. R.; Masle-Farquhar, E.; King, C.; Skvortsova, K.; Weatheritt, R. J.
Show abstract
The X chromosome carries an unusually high density of immune genes and is a major contributor to sex differences in immune function and autoimmune diseases. In females, X-chromosome inactivation (XCI) has two major functional consequences: it shapes X-linked gene dosage through XCI escape and determines the cellular exposure of heterozygous X-linked variants through XCI skew. Yet because XCI creates a mosaic of cells expressing different parental X chromosomes, these properties have remained largely inaccessible in individual women, becoming measurable only where XCI is non-random or after aggregation across large cohorts. Consequently, how X-linked variation contributes to sex-biased immunity and differs between individual women has remained unresolved. Here we present scDaisyChain, a graph-based framework that reconstructs chromosome-scale X haplotypes directly from heterozygous SNPs and single-cell long-read transcriptomes. scDaisyChain achieves near-ground-truth accuracy in highly polymorphic mouse hybrids and shows strong concordance with orthogonal long-read whole-genome phasing in human samples. Applied to peripheral blood immune cells from healthy women, it reveals a lineage-specific escape program in which lymphoid cells escape XCI more broadly than monocytes, with corresponding gains in the inactive X chromatin accessibility and female-biased expression. Lineage-specific skew further alters the proportion of cells expressing each heterozygous X-linked variant, a property we term variant exposure. Predicted deleterious variants are preferentially found in low-exposure states, exemplified by a splice-altering TLR8 variant expressed in few cytotoxic T cells. In rheumatoid arthritis (RA), the monocyte compartment - which has the lowest escape in health - shows reproducible inactive X dysregulation converging on a trained-immunity programme linked to disease flare and synovial macrophage activation, with elevated escape of IL13RA1 and HDAC8. These findings establish lineage-specific escape, skew and variant exposure as quantifiable, patient-resolved determinants of sex-biased immune gene dosage and X-linked variant penetrance in health and autoimmune disease, resolving a dimension of female biology that has been previously inaccessible in individual donors.
Dessart, M.; Luff, S.; Smith, L.; Sunman, H.; Vinauger, C.
Show abstract
Circadian clocks enable mosquitoes to anticipate recurring environmental variations and coordinate behaviors critical for survival and disease transmission, such as locomotion, reproduction, host-seeking, and blood-feeding, with times of day when performance is maximal. In Aedes aegypti, locomotor activity follows a robust diurnal rhythm shaped by endogenous circadian clocks and environmental cues, among which light has been shown to be the primary source of temporal information. While early studies established the role of light in regulating locomotor activity, behavior, oviposition and pupation, it remains unclear which features of a light cycle drive changes in circadian rhythms. This question is increasingly relevant as Ae. aegypti is frequently exposed to artificial and dynamic lighting conditions in urban environments. Here, we investigated how transient changes in light schedules influence circadian rhythms in locomotor activity by systematically manipulating the timing, duration, and direction of light exposure. Using a high-throughput assay, we tested over 1900 individuals, including wild-type and timeless knockout mutants, and showed that a single day of al tered lighting is sufficient to induce robust phase shifts, with no evidence of masking effects. A 6-hour light pulse was sufficient to re-entrain mosquitoes regardless of the timing of the pulse, and phase shifts were primarily driven by the offset time of the light pulse, indicating that light-offset acts as a major zeitgeber. Together, these findings challenge conventional assumptions about the timescale of circadian synchronization and highlight the remarkable plasticity of mosquito behavior in response to anthropogenic light. Eventually, these effects could explain the rapid adaptation of the species to urban environments and have potential consequences for disease transmission dynamics.
Yamaguchi, K.; Uchida, K.; Hiraiwa, M.; Fukano, Y.
Show abstract
Citizen science observations are abundant, but conservation requires turning uneven records into reliable predictions and directing new surveys to where information is missing. We developed a biodiversity platform for Japan that is updated monthly and integrates 2.32 million records to predict 8,297 species across seven taxonomic groups. Shared representation models outperformed species-specific models in four groups and extended predictions to species with few records. Five independent datasets, including structured monitoring, environmental DNA and complete forest inventories, confirmed that the models ranked observed species and occupied sites above alternatives, with median AUCs of 0.724 to 0.894 across sites and 0.650 to 0.841 across species. For any user-selected area, the platform returns candidate species, distribution predictions, a biodiversity map corrected for uneven observation effort, a conservation priority map for native species and a map recommending where to survey next. This map highlights places where species with few records are predicted to occur despite limited sampling. Independent observations showed that areas ranked highly by this predicted potential contained many such species, indicating that model predictions can help direct surveys toward knowledge gaps. New observations are incorporated into monthly updates, creating a national feedback system connecting citizen science, local conservation decisions and future surveys.